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ABSTRACT 


Phenological data as reported- at the crop reporting district (CRD) level, 
and environmental data from the National Climatological Center appear to be 
viable sources of information for use in development and testing of an ad- 
justable crop calendar model for winter wheat. These data were utilized for 
this study, in which generalized least squares techniques were applied for 
parameter estimation of functions to predict winter wheat phenological stage, 
with environmental values as independent variables. The independent variables 
investigated included daily maximum temperature ( , daily minimum tempera- 
ture (T ), daily daylength (D t ), and daily precipitation (P ). 

After parameter estimation, tests were conducted using independent data. 
From these tests, it may generally be concluded that exponential functions 
have little advantage over polynomials. Precipitation was not found to signi- 
ficantly affect the fits. The Robertson "triquadratic" form, in general use 
for spring wheat, was found to yield good results, particularly for the 
emerge-Joint interval, but special techniques and care are required for its 
use. In most instances, equations with nonlinear effects were found to 
yield erratic results when utilized with averaged daily environmental values 
as independent variables. 

The particular combination of gr.'wth rate functions recommended, within 
the various phenological intervals is: 

Plant-Emergei R = 0. 087870 - 0.000^5898 \? M 


ii 


Emerge-Joint : 


Joint-Head: 


Head-Soft Dough: 
Soft Dough-Ripe: 


R * 0.07X78'* (D l - 0.2796} * 

{(4.6869 x 107 l ‘)(T x + 2.3876) - 
(4.6618 x 10‘ 6 )(T x + 2. 3976) 2 + 

(2.3943 x 10"S(T x + 2.3876) + 

(1.7055 x 10" 5 )(T M + 2.3876) 2 ) 

R * exp {-6.22319 + 0.12008 T x 
«■ 1.0337 (D l ‘12) - 0.038591 T x (D l - 12)) 
R = 0.037980 ♦ 0.006ll45 T m 
R = exp {-17.6837 + 0.53654 T x 
+ 0.33366 T M + 2.6105 (D l - 12) 

-0.0080340 T x 2 - 0.011446 T^ 2 


-0.45441 (D l - 12) 2 }. 

Variance propagation studies indicate that the results obtained utilizing these 
functions are consistent with the statistical uncertainties associated with the 


data used for parameter estimation and testing. 
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INTRODUCTION 

The Large Area Crop Inventory Experiment (LACIE) utilizes winter wheat 
crop calendar, or biometeorological time scale, information in several ways. 
Among these, analyst interpreters use the information in training field 
selection and other researchers use the information in yield modeling efforts. 

To date the needed winter wheat crop calendar information as been pro- 
vided in one of three ways. In the first, historical records are used to com- 
pute an average date at which the crop will reach a specific stage, A 
"spread", computed statistically from the data, is included to account for 
yearly variations. In the second, functions found by Robertson (l) are 
modified for use in winter wheat. These modifications are necessary be- 
cause the Robertson's model was generated for spring wheat in the Canadian 
prairie provinces. The most commonly used procedure used to modify the 
Robertson's model has been to calculate "multipliers" for winter wheat (4), 

A third technique (5) utilizes the "normal" winter wheat crop calendar, 
modified for "location specific" effects, as indicated by long term tempera- 
ture averages . 

These techniques suffer from serious limitations. In the first, 
variations in planting dates, and the large variations inherent whe., aver- 
aging over many years, means that considerable uncertainties result in the 
estimated dates. For the second, it must be emphasized that Robertson's 

t 

model was generated for spring wheat, whose growth environment differs 
*Nwbers in parentheses throughout this report refer to the reference list in 


the back. 
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considerable from that of winter wheat. Spring wheat temperatures in General 
are on the rise throughout the growth and ripening of the crop. Daylengths 
in general increase during emergence and subsequently decrease. Winter wheat, 
on the other hand, is planted and emerges during a period of declining temper- 
atures and daylengths, undergoes a long winter dormancy period, and then 
completes its growth under conditions of increasing temperature and daylengths 
after spring greenup. The third technique yields good results, but does not 
yield a universal model applicable at locations worldwide, being dependent 
upon long term temperature records. 

For these reasons it was felt desirable to seek an acceptable crop 
calendar for winter wheat specifically, based upon phonological reports, ip 
which weather and environmental variables could be used to predict winter 
wheat growth stages. This was the overall goal of the research conducted. 

With this goal in mind the following research objectives were formulated': 

a. to define data sources which may be utilized in the formula- 
tion of such a model; 

b. to generate computer programs which may be used to investigate 
various mathematical model; 

c. to identify functional forms which are appropriate for des- 
cribing winter wheat growth stages; 

d. to obtain reliable estimates of parameters for these func— 

i * 

tions, and variance estimates of these parameters; 

e. to identify and investigate auxilliary problem areas, such as 
winter dormancy criteria; 

f. to test the parameter estimates obtained using independent 


data sets ; and 
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g. to perform variance propagation studies to assess expected 
variations in applying the -models, and compare these with 
variations realized during testing. 

In the first phase of the investigation (6), activity was centered upon 
gathering phenological and meteorological data at the Crop Reporting District 
(CRD) level, reduction and editing of i;M?» data, and application of generalized 
least squares techniques for parameter estimation. Evaluation of the resulting 
functions was limited to statistical testing of residue 1 ' tnd variances of 
residuals resulting from the least squares estimation procedure. 

During the second stage of the investigation, increased data were gathered 
and utilized in the parameter estimation programs, test programs were gen- 
erated and data produced for use in these, and variance propagation studies 
were conducted. 


FUNCTION ESTIMATION 


The objective of the modeling process was to obtain a relationship pre- 
dicting phenological stage numbers as a function of environmental variables. 
The phenological stages and the associated stage numbers used in this study 
are as follows: 


Stage 

Planting 

Emergence 

Jointing 

Heading 

So ft -Dough 


Number 

0 

1 

2 

3 

k 


Ripe 


5 


u 


The environmental variables selected were the same as those used by 
Robertson (l). These variables are daily maximum temperature, daily minimum 
temperature, and daylength, used to represent photoperiod. In addition, 
precipitation was investigated as a possible independent variable, in an 
attempt to assess the possible importance of moisture in such a model. 

Data Sources 

Phonological data at the Crop Repoi'ting District (CRD) level were used 

in this study. Depicted in Figure 1 are the crop reporting districts used 

in the United States. Historical data from some Great Plains, Midwest, and 

Rocky Mountain states were obtained from NASA's Johnson Space Center, Earth 

Observations Division. Data were gathered and included in the data sets for 

least squares parameter estimation. A summary of the number of location- 

years included for each stage is shown; 

Stage Location-Years 

Plant-Eraerge (O-l) 37 

Emerge- Joint (l-2) 53 

Joint-Head (2-3) 79 

Head-Soft Dough (3-4) 6 5 

Soft Dough-Ripe (4-5) 62 

For each location-year of data, meteorological data were averaged from several 
stations within the crop reporting district. Specific location-years used in 
the least squares program are summarized in Table 1. The dates used for 
modeling were the 505?, or median dates, converted to Julian Day Numbers, at 
which the crop reached each stage. Table 2 shows a typical report of the 
dates at which the crop reached each stage for a single year in one state. It 
should be noted that, throughout the United States, there exists a great 
variability as to winter wheat phenology reporting. As may be seen from 
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Table 1 

Summary of Location Years Used 
For Least Squares Parameter Estimation 


Staple 

State 

Crop Year 

C • R • D • 

0-1 

Colorado 

1972 

1,2,7 * 8,9 



1975 

1 , 2 , 6 , 8, 9 


Idaho 

1975 

1,9 


Oklahoma -3- 

1965 

2 ,3 ,*4 , 5 , 6 ,7 , 8 ,9 



1966 

1 , 2 , 3 , *4, 5 , 6 , f , 8,9 



1971 

1,2,3,14,5,6,7,8 

1-2 

Colorado 

1972 

1 , 2 , 6 , 7, 9 



1974 

2 , 6 , 8, 9 


Idaho 3 

1975 

1,9 


Montana *’ 2,3 

1971 

1.2,3 



1972 

1,2 


North Dakota *’ 2 

1970 

• 1,*4,7 ,8,9 



1971 

1,2,14,7,8,9 



197*4 

1 


Oklahoma 

1965 

2,3,14,5,6,7,8,9 

• 


1966 

1 , 2 , 3, 4, 5 , 6 , 7, 8 , 9 



1970 

1 , 2 , 3 , 4, 5 , 6 , 7, 8 

2-3 

Colorado 

1972 

2, 6, 7 , 9 



197*4 

2,6,9 


. Idaho 3 

1973 

1,9 



1975 

1,9 


Kansas 

196 k 

1 , 2 , 3, 4, 6 , 7, 8 , 9 



1965 

1,2,3, 5, 6 , 7, 8, 9 



1968 

1 , 2 , 3 , 4, 5 * 6 , 7 , 8,9 


Montana 3 

1971 

1,2 



1972 

1,2,3 


North Dakota 

1970 

1,4,7, 8,9 



1971 

1 , 2 , 4, 7, 8, 9 



197*4 

1 


Oklahoma 

1965 

1 , 2 , 3 , 4, 5 , 6 , 7 , 8,9 



1966 

1 , 2 , 3, 4, 5 , 6 , 7 , 8, 9 



1969 

1 , 2 , 3 , 4, 5 , 6 , 7, 8 
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Tabic 1 
(Continued) 


Stage 

State 

Crop Year 

C.R.D. 

3-4 

Colorado 

1974 

2,6,9 


Idaho 

1973 

1,9 



1975 

1,9 


Kansas 

1964 

1,2, 3,4, 6, 7 ,6,9 



1965 

1,2, 3, 4, 5, <5, 7, 8, 9 



1966 

1,2, 3, 4, 5, 6, 7, 8, 9 



1968 

1.2, 3, 4, 5, 6, 7, 8,9 


Missouri 

1970 

9 



1972 

9 


Montana 

1971 

1,2,3 



• 1972 

1,2 


Oklahoma 

1965 

1,2, 3, 4, 5, 6, 7, 8 



19 69 

1,2, 3, 4, 5, 6, 7, 8 

4-5 

Colorado 

1972 

6 


Idaho 

1973 

1,9 



1975 

1,9 


Kansas 

1964 

1,2, 3, 4, 6, 7, 8, 9 



1965 

1,2, 3, 4, 5, 6, 7 ,8,9 



1966 

1,2,3, 4, 5, 6, 7, 8, 9 



1968 

1,2, 3, 4, 5, 6, 7, 8, 9 



1969 

2, 3, 4, 5, 6, 7, 8, 9 



1974 

1,2, 3, 4, 5, 6, 7 ,8 


Missouri 

197o 

9 



1972 

9 


Montana 

1971 

2,3 



1972 

1,2 


*Data included only for winter dormancy through spring greenup. 

1. ) "Acceptable Stand" reported. Emergence estimated by interpolation, 

2. ) Emergence not reported. 

3. ) "Boot" reported. Jointing estimated by interpolation. 


• < 


Table 2: Typical Report of Winter Wheat Phenology for CRD's 


.l:__ . KANSAS 
C> op Yejr: _ 72, _ 


GROWTH STAG'i DUTIES reft 50 J DEVELOP. ’EM 
Month-Day -Year 


U 

CfV‘_\ C.i 1 ?* 

Eprin j ; 


as. 

1 ! 

car 

,'t:. 

Pl-lliti.ip 

Dr. to y 

Erccrgancs 
Dale y 

.Jointing 

Zf 

Hc.’C;rq 
Dote 2/ 

10 N.W. 

m 

Sept. 20 

* 

May 6 

May 24 

» r 
o 2 

20 W.C. 

Sept. 20 


May 2 

May 19 

c > 
>g 

30 S.W. 

Sept. 21 


April 20 

May 9 

-?o> 

40 N.C. 

Sept. 29 

H 

•§ 

H 

May 1 

May 20 ‘ 


30 C. 

Oct. 2 

«3 

< 

April 26 

May 14 


60 S.C. 

Sept; 28 

: S 

: z 

April 15 

May 6 


70 N.E. 

Sept. 29 

• 

• 

• 

• 

April 12 

May 22 


80 E.C. 

Oct. 1 

• 

• 

• 

April 10 

May 16 


90 S.E. 

Oct. 8 

• 

• 

• 

m 

April 15 

April 20 


Sofu 
Coup.; 
Dara 3/ 


June 18 
June 12 
June 5 
June 11 
June 6 
June 1 
June 11 
J"Tii } 
May 24 


1/ 0c cit .*1»h 50% cr crup ii. th* CPJ- was pTvi.rcl or or:?rr.c1, ; espset vc*»y. 


2- •' *' 
J, r *;M\ 


' c r. o/ C :v.- >. r-~ : C’3 V'! bs-»*»«i to }•??.': o.* •*:»*}, * ;='~t ?iy. 


Rip a 


Late 4/ 


June 30 
June 30 
June 21 
June 27 
June 21 
June 14 
June 29 
June 24 
June 13 


July 1' 
July - 
June 2.- 
July : 
June 2- 
June 2c 
June 30 
June 29 
June 2! 


i I 


r jV c U ' ECT UfT. tit Wt-ST T-Cfl 'bwCJ * {',Ur.-l.C 40*. ‘i* Cl 


- x ! K • o7 '.re.. ir is In vib - t.’S.V -ru^h) - 

‘...C-;** . •* i .« i- ;.p. *«<;«•■> c- . 


6t *,u i c* r-' r.f erno ir, ch~ C'5 fr .ir. fe.-ves* •*& either as stmOi\q c^air or t-t of <w4.*. 
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Table 1, few states report all the phonological stages, and several report 
stages which differ from the standard stages, For example, Oklahoma reports 
"acceptable stand" rather than emergence, and both Idaho and Montana report 
"boot" rather than Jointing. 

As independent variables, values of daily maximum and minimum temperature, 
daily precipitation, and daylength may be used. Values for temperatures and 
precipitation are available on a daily basis from the National Climatological 
Center, NOAA, Asheville, North Carolian. Records were obtained for the states 
and years for which phenological data were available. Tables 3 and k show 
typical listings of temperatures and precipitation, respectively, as reported 
by the National Climatological Center, The value of daylength may be cal- 
culated as a function of latitude and Julian Day Number. The empirical equa- 
tion obtained by Stuff (S) was used for this purpose in this investigation in 
which a tangent function is utilized for latitudes less than 40°, and a 
tangent-squared function for latitudes above Uo°. 


Functional Forms 

These data were used in general least squares procedures to estimate 
parameters for functions of the form 


R 


f( 


H X , 
p,l’q,r 


Eq, 1 


in which R represents daily rvowth rate, H represents a vector of parameters, 
variable with the model chosen, estimated through the generalized least squares 


L J 
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algorithm, and. X represents the vector of independent variables chosen, in this 
case, 

ft = [ T x T M D L P r] t 

where, T Y = daily maximum temperature in °C. , 

T„ = daily minimum temperature in °C. , 

M 

D. = day length (representing photoperiod) in hours, 

L 

* daily precipitation in millimeters. 

The specific functional forms selected were of three classes. The first 
was a polynomial of the form 

,=I 1* H 2 T X + l, 3 T M + VV 12 > + Vr 
+ V*® + H 7 T M 2 + H 8 ( V 12)2 + V/ 

Eq. 2 

+ H 10 T X^ D L" 12 ^ 4 ir il T M^ D L -12 ^ + II 12 P r" D ir 12 ^ 

* Wm + Wr + H l 5 Vr- 

The second functional form investigated was the exponential 
R = exp (H x + H 2 T x + H ? T m + H^-12) + 

♦ H 6 T X 2 + H 7 T M 2 + H 8 (D 1- 12)2 + V/ 

Eq. 3 

+ H 10 T X ( V 12) + H 11 T M ( V 12) + H 12 P r ( V 12) 


+ H 13 T X T M + H ll» T X P r + H 15 T mV' 
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For these functions, several cases were run using computer programs 
generated during the first phase of the contract, in which various subsets of 
the equations were chosen. In all cases, terms containing daylength (D^) 
were omitted for stage 0-1. 

The third functional form investigated was the "Robertson model" (l), of 
the form 

R = {H 2 (D l - H 1 ) + H 3 (D l - H^ 2 } • {H 5 (T x - iy 

Eq, 1+ 

+ a 6 (t X - H U )2 + H 7 (t M ‘ V * h B (t m - V 2) 


for all stages except the plant-emerge interval. For this stage, a function 
of the form 


R = {Hg(P r - H 1 ) + H 3 (P r - H 1 ) 2 } • {H 5 (T x - H^) 

+ V T X ' \ )2 + H 7 (T M - V + H 8 (T M " V 2) 


Eq. 5 


was chosen. 


Least Squares Techniques 

For this investigation parameters for the above functions were estimated 
using the general constrained minimum least squares approach as published 
by Mikhail (3). The method, called "Simultaneous adjustment of observations 
and parameters," allows maximum flexibility for investigations of alternate 
functional forms and provides an extemely powerful tool for determining 
parameters of these functions and estimating variances of the parameters 


lU 

obtained. The following discussion provides a brief summary of the technique. 


For a more detailed discussion, the reader is referred to Mikhail (3). 
Given a vector of n original observations 


= [ X l’ X 2 ,,,,X n] 


n,i 


Eq. 6 


and a set of condition equations of the form 


r,l u,l n,l r,l 


Eq. 7 


In which there are r equations, X is the vector of n adjusted observations, 

H is the vector of u parameters to be estimated, and is the null vector, 
then the following may be stated. 

Since, in general, these equations will be nonlinear, a Taylor’s series 
approximation may be written as 


A V _B A_ _ F° 

r,n n,l r,u u,l r,l , Eq. 8 


in which 




B 


3F 

3H 


H° 



Eq. 9 


are the Jacobian matrices of the functions with respect to the observations 
and parameters, respectively, and evaluated at observed values and parameter 
estimates. Further, 



Eq. 10 


i ft* 
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Is the evaluation of the functions with observed values and parameter 
estimates, V is the vector of residuals to be applied to the observed 
quantities, and A is the vector of corrections to be applied to the 
parameters. 

The solution proceeds using the least squares condition in which 


V t P V 
l,n n,n n,l 


Eq. 11 


is to be minimized. P is the weight matrix of the original observed data. 
If a priori estimates of variance are known, then 



n,n n,n 


Eq. 12 


in which is the a priori variance/covariance matrix of the observed 
quantities. If this variance/covariance matrix is unknown, then P represents 
some assumed weight matrix, and may be written as 


P_ = Eq. lo 

in which Q t q represents some assumed cofactor matrix of relative variances 
and covariances. Since, in general, the A Jacobian matrix in Equation 8 will 
not be square, and hence, the linearized condition equations cannot be solved 
directly for the residuals and substituted into the least squares condition 
of Equation 11, the more general method of constrained minimum is used. In 
this method, a vector of r Lagrongian multipliers is defined. 
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h 




• • 



Eq. 14 


Then, the following scalar is minimized rather than 11. 

♦' = / PT- 2jv£(A V + B A - F°) Eq. 15 

Minimizing by taking partials of thi3 scalar with respect ot the observations 
and parameters and augmenting with the original linearized condition equa- 
tions resulting in the following total set of normal equations: 


P A t ♦ 


V 


'_4_' 

n,n n,r n,u 


n,l 


n ,1 

A 4> B 


A 


F° 

r,n r,r r,u 


r,i 


r ,1 

4> B 4> 


A 



u,n u,r u,u _ 


_ u,l ^ 


U,1 


Eq. l6 


This set may be solyed by partiti'r ling with the following results 


A _ _T_ 

U,1 ” U,U U,1 


Eq. 17 


-1 


K L = (A£ L oA b ) (-B A + F ) 


'Eq. 18 


V = fi L o h. K x 


Eq. 19 


r 
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in which, 

W m P'iVi*’' 18 ] 

^=b\a V a‘)-V. 

These equations are solved in an iterative manner until convergence. 

After convergence, the reference variance may be computed by 

2 I* El 
S o ~ r-u 

Variance propagation techniques may be used to show that (3), 

S = s 2 N" 1 Eq. 23 

“X o — 

yields an estimate of the a posteriori variance/covariance matrix for the 
parameters . 

The reference variance s q represents a measure of the "goodness" of the 
fit, and may be used in F tests to determine the statistical significance in 
adding or deleting terms from the functional form under investigation* or 
in comparing different functional forms. 

For the problem at hand the function of Equation 1, formulated in terms 
of Equation 7» would be for each data point 


Eq. 21 


Eq. 22 


Eq. 20 


F = R - f(H,X) = 0 


Eq. 2k 
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In which r * 1, n * 5* and u would depend upon the number of parameters 
required for the functional form selected. For the case using the general 
least squares procedure, 

sfl ' [ R T X t m :l L P rJ E? - 25 

uTl * [ h l’ h 2 ,,,,h u] t - Eq * 26 

IMPLEMENTATION 

Efforts were made to formulate and execute a data reduction system sub- 
ject to the constraints of the data forms available and consistent with the 
data requirements for application of the least squares techniques discussed in 
the previous section. Computer programs were written for the polynomial and 
exponential function forms using both the variation of parameters and the more 
general least squares techniques. Special techniques were required in treating 
the Robertson model. 


Data Reduction and Prenar ation 

Although the objective of the research was to model daily growth rates, R, 
as functions of environmental variables, no data could be obtained for growth 
rate on a daily basis. As may be seen from Table 2, phenological data re- 
ported consisted at best of dates at which the crop reached each phenologi- 
cal stage. With this data, only an average growth rate could be computed. 


of the form of 


in which n = the number of days in the i-th phenological interval for the 
location/year. This average growth rate represents an average, not only 
over the time period for the interval, but also a spatial average for the crop 
over the geographic region representing the crop reporting district. 

Environmental data were available on a daily basis. However, it would be 
inconsistent to apply daily environmental data in the least squares modeling 
process when only averaged growth rates were available. Therefore, in order 
to supply environmental data consistent with the growth rates available, the 
following pre-processing steps were performed. 

First, within each crop reporting district several meteorological re- 
porting stations were selected. Up to five stations were used. Temperature 
and precipitation data were recorded for each station on a daily basis through- 
out the reported phenological stage interval from the climatological records. 
Table 5 represents a typical data reduction form used for this purpose. 

This information was keypunched onto cards and input into a data re- 
duction program. This program computes the average position of the crop 
reporting district (latitude, longitude, elevation) and the associated 
statistics, converts the meteorological data from English to SI units, 
averages meteorological values (T^,T^,P r ) for the location/year, and computes 
the associated variances. A listing of this data reduction program with 
output from a typical run is shown in (6), For all of the phenological inter- 
vals, with the exception of the emerge-Joint interval, the results for each 
location/year would be a single data point for input into the least squares 
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program, containing averaged growth rate, maximum temperature, minimum temper- 
ature, daylength, precipitation, and positional information for that location/ 
year. 


Winter Dormancy and Spring Greenup 
Special considerations were necessary for the emergy-Joint interval. 
Because this interval stretches over the long winter dormancy period, no 
single averaged growth rate value could he considered to reasonably repre- 
sent the entire interval. For the growth rates the interval was broken into 
three sections, as depicted in Figure 2. A fall emerge-dormant period must 
be defined. The crop is then assumed to have completed one half of its inter- 
val growth in this period, and a growth rate of 


R E 


0.5 

N E 


Eq. 28 


was assigned for this period. Over the winter dormancy period, N^, a growth 
rate of zero was assumed. Lastly, a spring greenup-J oint interval is defined, 
and a second average growth rate computed as 



Eq. 29 


These rates, and the environmental data averaged on a monthly basis were used 
to define several data points for each location year in this interval, 

Therefore, it was found necessary to define some criteria for winter and 
spring. For this purpose, a temperature criterion was selected in an attempt 



Figure 2: Emerge- Joint Growth Rate Assumptions 


1 
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to predict the onset of winter and spring, Growth rates and temperatures used 
for stage 1-2 in the least squares programs were input to the Cat e-Liebig (7) 
critical level analysis program. Figures 3 and represent plots of growth 
rate versus temperature for the winter and spring runs, respectively. Criti- 
cal levels computed by the algorithm were h.69°C. (l40.l*°F.) for the onset of 
winter dormancy, and 6.lk°C. (43.1°F.) for subsequent spring greenup. 

Therefore, for winter, the criterion that 10 days of average temperature 
below 4,69°C. was used. For spring greenup, the criterion of 10 days above 
6.lh°C. was used. 


Computer Runs 

The functional forms of Equations 2 and 3 were programmed for parameter 
estimation using the generalized least squares algorithms, A listing of this 
program, with typical output is shown in (6). For the investigation, func- 
tional forms were selected and computer runs made using selected terms from 
Equations 2 and 3. For the polynomial, runs were made using only the constant 
term, using single variable linear equations in D^, T^, and T^, using linear 
combinations of these variables with and without precipitation, and various 
selected combinations of linear with interaction terns and linear with 
squared terms. The same set of combinations was run for the exponential 
function. 

Specific features of the programs which may be of interest include the 

i 

fact that any functional form may be accomodated by changing a single statement 
in the EF function subprogram, Linearization is accomplished by numerical 
evaluation of partial derivatives using forward differences. Even greater 


i 
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Figure 3 

Plot of Data Input to Cate-Liebig 
for Winter Dormancy Criteria 


»•&* 06 TA ►Ok VfckttnLc 






Plot of Data Input to Cate-LieMg 
for Spring Greenup Criteria 
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flexibility is produced by the introduction of a parameter selection index, 
IPAR, which allows the inclusion or deletion of parameters within the func- 
tional form being investigated. 

The use of such numerical techniques, which allows great flexibility, 
is not without danger. It was found in some instances that convergence of the 
solution did not occur when it was to be expected due to the approximations 
inherent in these numerical techniques. These problems occurred only for the 
highly nonlinear functions and is not felt to be a serious problem for reasons 
to be discussed subsequently. 


The Robertson Model 

As discussed previously, the "triquadratic" model advocated by Robertson 
has much to recommend it. The model. Equation 4 , can accomodate both non- 
linear effects and interactions with a minimum number of parameters. However, 
the model requires special considerations in estimating parameters. As 
pointed out by Robertson in his paper (l), the parameters for the model are 
not independent. An attempt to simultaneously estimate all of these para- 
meters using the generalized least squares technique verified this, for any 
attempt at such modeling results in a near-singular normal equation matrix. 

Robertson utilized a parameter transformation technique. By expansion 
of each multiplicative factor separately, polynomials result which have the 
same number of coefficients as original parameters within each factor. 
Estimation may then be carried out by holding the coefficients of the trans- 
formed daylength parameters, for example, and using regression to estimate 
transformed temperature parameters. These temperature parameters may then 
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be held, and transformed daylength parameters found by regression. This 
successive substitution technique is cycled until no change in the parameters 
result, and the parameters are transformed back for the original equation 
form. The important point here is that the temperature parameters are 
estimated independently from those for daylength. 

The same technique was implemented more simply using the generalized 
least squares techniques discussed previously. The procedure is identical, 
except no transformation is required, since the method is not limited to 
regression. The method proceeds as follows. First, parameter approximations 
are assumed for the daylength factor. These are held as constants, and the 
general least squares algorithms are applied in an iterative manner until 
convergence for the temperature parameters. These are then held, and 
another computer run made to estimate new values for the daylength para- 
meters, using these generalized least squares algorithms. The process is 
repeated in successive computer runs until no change is noted in the 
parameters. 

• Parameter Estimates 

Parameters generated from the least squares programs are summarized in 
Tables 6 through 16. These parameters were then used as input for the sub- 
sequent testing and variance propagation programs. An important feature of the 
generalized least squares modeling technique is that a-posteriori estimates of 
the variance-covariance matrix for the parameters result from its application. 
These may then be used in subsequent variance propagation studies. 
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Table 7 

Polynomial Coefficients Resulting from Least 
Squares Fit, Stage 1-2, Based on 53 Location Years 



Const . Only 

Lin Only 

Lin T fJ Only 

Lin D Only 

ii 

H (X10~ 3 ) 

1*. 859397 

0. 8615316 

5.7987^6 

5.634711 

H 2 (X 10 _li ) 

0 . 

3.4423 

0 . 

0 . 

H.(XlCf U ) 

_li 

H U (X10 ) 

0 . 

0 . 

0 . 

0 . 

3. 91*45 
0 . 

0. 

9.9291 

h 5 (xio' u ) 

0. 

0. 

0. 

0. 


T x’ d l 4 t x°l 

V d l & Vl 

Lin & Sq v/o P^ 

Lin & Sq v/ P f 

IL^XIO -3 ) 

HgCXlO" 14 ) 

2.413333 

2.6036 

6.192233 

0. 

4.749293 

1.3716 

3.635508 

1.9310 

H 3 (X10 _lt ) 

0. 

2.8102 

2.1683 

1-6575 

H^XIO'S 

7.6379 

0.78146 

1.3601 

1.1765 

h 5 (xio _ 

0. 

0. 

0. 

5.7270 

H g (X10 -5 ) 

0. 

0. 

0.10771 

4.7487 

H (X10 - 5 ) 
HgCXIO -5 ) 

0. 

0. 

0. 

0. 

0.83089 

-45-599 

0.93968 

-47.176 

H (X10~ 5 ) 
H 10 (no-5) 

0. 

- 5.1193 

0. 

0. 

0. 

0. 

-7-7746 

0. 

H u Oao-=) 

0. 

-6.3540 

0. 

0. 

h 12 (xio- ) 

0. 

0. 

0. 

0. 

h 13 (:qo- ) 

0. 

0. 

0. 

0. 

H^tXlO" 5 ) 

0. 

0. 

0. 

0. 

h 15 (xio -5 ) 

0. 

0. 

0. 

0. 


Lin P r Only 

Lin Ty., T m , D l Lin 

V V D L’ P r 


3-958185 

2.755673 

2.222651 


0 . 

2.1425 

2.3935 


■ 0 . 

1.5561 

1.2110 


0 . 

0.14262 

0.082501 


7-3297 

0 . 

1.2524 


Lin & Int v/o P 

Lin & Int v/ P 

r 

All v/o P r 

All v/ P r 

1.906250 

-3.485117 

14.44650 

-3-501598 

2.6840 

6.1158 

-13-507 

6.3881 

- 2.0061 

-5-5208 

18.229 

-4.0224 

6.1203 

- 9.6170 

-3.5812 

-25-772 

0. 

28.800 

0. 

47.081 

0. 

0. 

5.4l6l 

0.11238 

0. 

0. 

7-5573 

0.58963 

0. 

0. 

-35.466 

-44.423 

0. 

0. 

0. 

-13-640 

-3.6078 

5.9229 

3.2799 

16.239 

-9.5090 

-17-6415 

-7.2848 

-19.051 

0. 

21.548 

0 . 

26.190 

1.6780 

1.4762 

- 11.831 

0.36190 

0. 

- 20.016 

0. 

-28.796 

0. 

23.010 

0. 

31.537 


Table 8 

Polynomial Coefficients Resulting from Least 
Squares Fit, Stage 2-3, Based Upon 79 Location Years 



Const. Only 

Lin Only 

Lin Ty Only 

Lin D t Only 

Jj 

H 1 (Y10 -2 ) 

U . 5896 U 8 

2.782055 

2.966387 

3.519618 

h 2 (xio" 3 ) 

0 . 

0.74951 

0 . 

0 . 

H 3 (:acf 3 

0 . 

0 . 

1.8460 

0 . 

H, (XIO -3 ) 

0 . 

0 . 

0 . 

4.7037 

H (X10 ) 

0 . 

0 . 

0 . 

0 . 


T X’ d l 4 t x°l 

V d l & T A 

Lin & Sqw/o 

Lin & Sq v/ P 

H 1 (X10 -2 ) 

- 21.79267 

-2.129727 

53-03508 

42.06301 

h 2 (xio~ 3 ) 

10.621 

0 . 

-43.989 

-37-715 

h 3 (xio -3 ) 

0 . 

5.6399 

4.8274 

-0.77791 

h 4 (xio -3 ) 

159 . 42 

18.087 

26.863 

40.810 

H C (X10~ 3 ) 

5 

MX10 ) 

6 -4 

H„(X10 ) 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

8.4855 

0.30465 

3.6933 
7.9387 
1 . 4oi6 

H a (X10 - ^) 

H (XIO" ) 
H 10 (7- 10 -S 

0 . 

0 . 

-65.800 

0 . 

0 . 

0 . 

-45.496 

0 . 

0 . 

- 69.655 

4.3304 

0 . 

H (XK) ) 

J-i- L 

0 . 

-12.542 

0 . 

0 . 

H 12 ('X10- ) 

0 . 

0 . 

0 . 

0 . 

H..(X10" ) 
- L ^ -li 

0 . 

0 . 

0 . 

0 . 

H llt CX10 ) 

0 . 

0 . 

0 . 

0 . 

H 15 (X10 4 ) 

0 . 

0 . 

0 . 

0 . 


Lin P Only 
r 

Lin T x , T m , D l 

Lxn T x , T m , D l , P r 


3.906240 

4.519374 

1-723743 


0 . 

-1.4123 

- 0.061902 


0 . 

2.6964 

1.4967 


0 . 

4.8582 

5-0642 


2.7747 

0 . 

2.2210 


Lin & Int v/o P 

r 

Lin & Int v/ P 
r 

All v/o P 

r 

All v/ P r 

3.805418 

-21.47510 

- 308.7502 

66.21290 

0.80023 

-0.13505 

246.80 

-57-263 

- 15.437 

127-00 

-0.62885 

43-384 

62.333 

-493-18 

191.53 

-67.352 

0 . 

-8.782 

0 . 

-66.358 

0 . 

0 . 

-31-732 

11.624 

0 . 

0 . 

50.161 

11.535 

0 . 

0 . 

196.53 

-125.48 

0 . 

0 . 

0 . 

17-122 

- 36.190 

226.71 

-209.48 

63.259 

29.910 

-85-739. 

231-95 

-67-264 

0 . 

82.147 

0 . 

191.59 

5.4502 

-39-521 

-57-798 

-18.579 

0 . 

31.192 

0 . 

17-402 

0 . 

-77.920 

0 . 

-17-030 


Table .9 

Polynomial Coefficients Resulting from Least 
Squares Fit, Stage 3-4, Based Upon 65 Location Years 


• 

Const . Only 

Lin T^ Only 

Lin Tjj Only 

Lin D Only 

ii 

H 1 (X10' 2 ) 

4. 553196 

1*. 259662 

3.798018 

6.377203 

h 2 (xio -3 ) 

0 . 

1.0945 

0 . 

0 . 

H-Cxio -3 ) 

0 . 

0 . 

0.6lll*5 

0 . 

Hj^XlO -3 ) 

0 . 

0 . 

0 . 

-7.06i4 

H 5 (X10 -3 ) 

0 . 

0 . 

0 . 

0 . 


T x’ °L & t x d l 

T M’ d l 4 t m°l 

Lin Sc Sq w/o P r 

Lin Sc Sq v/ P 

H-jUlO -2 ) 

-7.890370 

-1.21*0937 

10.15888 

10.43590 

Hp(X10~ 3 ) 

5.1*959 

0 . 

6.71*31+ 

5.8296 

H 3 (X10“ 3 ) 

0 . 

6.9156 

-O .67838 

-1.2164 

H, (X10~ 3 ) 

1*3.655 

21.265 

110.78 

- 10.586 

H 5 (X10" 3 ) 

0 . 

0 . 

0 . 

0.41741 

H/-(X10~^J 
0 » 

II T (X10~ 4 ) 

HgCXlO" 1 ') 

000 

odd 

- 1.0120 

0 . 11*672 

186.18 

-0.79778 

0.33175 

179.47 

H (X10~^) 

0 . 

0 . 

0 . 

0.21120 

9 - 1 * 
H io (xl0 . } 

-19.733 

0 . 

O. 

0 . 

H.-CXIO" 4 ) 

0 . 

- 26.766 

0 . 

0 . 

H 12 (X1 °" ,) 

0 . 

0 . 

0 . 

0 . 

H_(X10 ) 

0 . 

0 . 

0 . 

0 . 

H lU (X10 ) 
h 15 (xio ') 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 


Lin P r Only 

Lin T x , T m , 

Lin V V V P r 


4.181928 

6.508311 

5.278253 


0 . 

0.0054597 

0.34962 


0 . 

- 0.072310 

-0.30235 


0 . 

-7.2799 

-6.0371 


1.0661 

0 . 

0.77576 


in Sc Int v/o P 

r 

Lin Sc Int v/ P^ 

All v/o P 

r 

All v/ P r 

- 3 O.OOI 6 O 

25.8450 9 

- 8 . 343650 

- 26.11011 

10.388 

9-2824 

14.557 

20.432 

23.351 

15.191 

12.271 

13.253 

60.871 

68.840 

-106.31 

-43-059 

0 . 

-4.1958 

0 . 

3.1980 

0 . 

0 . 

-1.1747 

-0.88633 

0 . 

0 . 

1.7558 

2.1305 

0 . 

0 . 

180.15 

142.28 

0 . 

0 . 

0 . 

-2.6497 

-12.119 

-17-243 

1.0105 

-17.363 

• -35.160 

- 21.568 

- 3.0211 

12.405 

0 . 

-30.425 

0 . 

-40.064 

-5-6042 

-3-439 

- 6.1099 

- 8.5883 

0 . 

5.169 

0 . 

2.5037 

0 . 

-0.86249 

0 . 

2.2688 


Table 10 

Polynomial Coefficients Resulting from Least 
ires Fit, Stage U-5, Based Upon 62 Location Years 
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Table 11 

Exponential Coefficients Resulting from Least 
Squares Fit, Stage 0-1, Based Upon 37 Location Years 


H 1 

H 2 (X10~ 2 ) 

h 3 (xio~ 2 ) 

Hj * _2 
h 5 (xio *) 

Const. Only 

Lin T^ Only 

Lin Only 

Lin D Only 

Li 

Lin P r Only 

Lin T x , T m , D l 

T x , V Bj,, P t 

All v/ P r 

-2.473168 

0. 

0. 

0. 

0. 

V d l * t x°l 

- 3 . 53661 U 

4.1*615 

0. 

0. 

0. 

T M’ d l * Vt 

-2.430835 

0. 

-0.55089 

0. 

0. 

Lin !c Sq v/o P y 

-2.473168 

0. 

0. 

0. 

0. 

Lin & Sq v/ P y 

-2.391737 

0. 

0. 

0. 

-5.3781 

Lin & Int v/o P- 
r 

-4.311131 

9.2731 

-4.8492 

0. 

0. 

Lin & Int v/ P 

T 

-4.536354 

10.315 

-5.5545 

0. 

2.0610 

All v/o P r 

H 1 

-3. 536611* 

- 2.430835 

-2.185849 

-2.048876 

-4. 318064 

-4.853732 

1.19859 1 * 

5.135680 

H 2 (X10 -2 ) 

1*. 4615 

0. 

-9.6644 

-12.374 

9.3044 

11.694 

-44.125 

-76.589 

H,(X10~ 2 ) 

0. 

- 0.55089 

-2.4199 

-2.9238 

-4.7443 

-10.534 

13.153 

14.582 

4 

0 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

H 5 (X10 -2 ) 

0. 

0. 

0. 

3.3268 

0. 

32.338 

0. 

- 19.235 

H 6 (X10~ 3 ) 

0. 

0. 

4.0728 

5.0017 

0. 

0. 

12.659 

19-093 

H 7 (X10 -3 ) 

0. 

0. 

-1.8834 

-2.3159 

0. 

0. 

-O. 78 B 02 

0.19905 

H g (X10~ 3 ) 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

H g (X10" 3 ) 

0. 

0. 

0. 

r 0.66833 

0. 

0. 

0. 

5.4270 

H. _(X10~ 2 ) 
10 „ 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

H-^XIO"*) 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

h 12 (xio~ 2 ) 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

H (X10 -2 ) 

0. 

0. 

0. 

0. 

-0.0044585 

0.17681 

-0.73298 

-0.78505 

h i4 (xio -2 ) 

0. 

0. 

0. 

0. 

0 . 

-1.1148 

0. 

1.1572 

h 15 (xio -2 ) 

0. 

0. 

0. 

0. 

0. 

-0 ,24084 

0. 

-1.0840 


1 

h 2 (xio -2 ) 

h 3 (xio" 2 ) 

H^XIO -2 ) 

h 5 (xio -2 ) 


1 

h 2 (xio‘ 2 ) 

H 3 (X10 -2 ) 

H^(X10~ 2 ) 

H 5 (X10 -2 ) 

H,(X10~ 3 ) 

6 

H (XIO J ) 
Hg(X10' 3 ) 
H g (X10~ 3 ) 
H 10 (X10 -2 ) 
H (XIO -2 ) 

h 12 (x 10 " 2 ) 

H 13 (X1 °" > 
K^CXIO- ) 

H ig (XlO _<i ) 


Table 12 

Exponential Coefficients Resulting from Least 
Squares Fit, Stage 1-2, Based Upon 53 Location Years 


Const. Only 

Lin T x Only 

Lin T„ Only 
M 

Lin D Only 
L 

Lin P Only 
r 

Un V V D L 

Lin T x , T m , D l , P r 

All v/ P r 

-5-326841 

0. 

0. 

0. 

0. 

c x’ °L 1 t x°l 

-6.169140 

6.3917 

0. 

0. 

0. 

V d l * t m°l 

-5.258371 

0. 

7.1782 

0. 

0. 

Lin & Sq v/o P y 

-5.241116 

0. 

0. 

12.361 

0. 

Lin & Sq v/ P r 

-5-483952 

0. 

0. 

0. 

11-799 

Lin & Int v/o P 

r 

-5-852482 

3.9636 

3.4320 

-4.4626 

0. 

Lin & Int v/ P 

r 

-5.963827 
4.5805 
2.6651 
-4.3100 
1.7655 
All v/o P 

T 

-6.352690 

-5.095987 

-9.589213 

-9.757012 

-5.994067 

-5-948517 

-6.151259 

-7.765557 

8.7016 

0. 

55-Ul*7 

57.556 

7.2653 

8.7771 

10. 552 

30-792 

0. 

7.0304 

2.3567 

2.7682 

29.750 

26,597 

40-575 

14.573 

68.805 

19.353 

10.293 

9-4373 

77.991 

73-821 

69.655 

33.945 

0. 

0. 

0. 

2.2833 

0. 

11.701 

0. 

49.615 

0. 

0. 

-15.036 

-15.699 

0. 

0. 

-0.69791 

-6.0984 

0. 

0. 

0.096239 

0.12464 

0. 

0. 

7.0172 

0-98346 

0. 

0. 

-74.122 

-71.564 

0. 

0. 

-79.068 

-71-219 

0. 

0. 

0. 

-7.9763 

0. 

0. 

0. 

-22.077 

-3.9052 

0. 

0. 

0. 

-4.7623 

-4.1519 

-3-9830 

-1.4824 

0. 

-5.6637 

0. 

0. 

3.6255 

2.1082 

3.6022 

-0.57119 

0. 

0. 

0. 

0. 

0. 

-O.I8591 

0. 

3-3078 

0. 

0. 

0. 

0. 

-1.4949 

-1.6837 

-2.3046 

-1.1756 

0. 

0. 

0. 

0. 

0. 

-2.4685 

0. 

-4.1073 

0. 

0. 

0. 

0. 

0. 

6.1144 

0. 

7.7660 
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Table 1*4 

Exponential Coefficients Resulting from Least 
Squares Fit, Stage 3-*4, Based Upon 65 Location Years 



Const. Only 

Lin T^ Only 

Lin T Only 
M 

Lin D_ Only 
I* 

H 1 p 

-3-0893*41 

-3.151508 

-3-2632*40 

-2.600222 

h 2 (xio" 2 ) 

0. 

0.123177 

0. 

0. 

h 3 (xio _2 ) 

0. 

0. 

1.14020 

0. 

Hj^XlCf 2 ) 

0. 

0. 

0. 

-19.117 

k 5 (xio~ 2 ) 

0. 

0. 

0. 

0. 


T X’ d l & t x d l 

t m’ d l t m d l 

Lin & Sq v/o P 
r 

Lin & Sq v/ P 

r 

H 1 

H 2 (X10' 2 ) 

-5-589585 
11. 369 

-*4.259738 

0. 

-1. 99*4119 
1*4.311 

-2. 03*4633 

13-55*4 

h 3 (xio -2 } 

? 4 (X10“ 2 ) 

0. 

88.113 

II4.I467 

*43.573 

-1.7818 

-2.2915 

-3.1076 

-2.2365 

H 5 (X10“ 2 ) 

H 6 (X10 -3 ) 

0. 

0. 

0. 

0. 

0. 

-2.1291 

0.98905 

1.8677 

h 7 (xio~ 3 ) 

0. 

0. 

0. *46881 

0.92379 

Hg(X10 -3 ) 

0. 

0. 

385.78 

380.38 

lyxio -3 ) 

H 10 (X10-=) 

0. 

- U. 1110 

0. 

0. 

0. 

0. 

0*27780 

0. 

a u (xio- ) 

0. 

-5.675*4 

0. 

0. 

h 12 (xio- ) 

0. 

0. 

0. 

0. 

H, 3 (X10 2 ) 
H lU (X10-=) 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

h i5 (xio d ) 

0. 

0. 

0. 

0. 


Lin P r Only 

Lin T x , T h , D l 

Lin T x , T h , P L , P 

r 

-3-168*43*4 

0 . 

0 - 

0. 

2.2172 

Lin & Int v/o P^ 

-2.612179 

0.22383 

-0.2391*4 

-19-837 

0 . 

Lin & Int v/ P^ 

T 

-2.923235 
0.89066 
-0. 677*49 
-1*4.739 
1.5822 

All v/o P r 

All v/ P r 

-11.2019*4 

- 10.07295 

-6.8*41593 

-3-03*4168 

2*4.875 

21.360 

32.728 

-27-975 

53-350 

3*4- *491 

31-721 

33.562 

1*42.98 

155-65 

-195-3*4 

-29. *453 

0. 

-6.7990 

0 . 

13.256 

0. 

0 . 

-2.0501 

15.376 

0. 

0 . 

* 4 . 1610 

3-6577 

0 . 

0 . 

355 . 9*4 

26*4.71 

0. 

0 . 

0 . 

-*4.3630 

- 3.0201 

- 3.8988 

-0.1975*4 

-*4.9320 

-7-69*40 

-*4.9288 

-0.99*436 

0.89*450 

0. 

-6.18*40 

0 . 

-5-9877 

-1-3176 

- 0 . 865 * 4 ** 

-1.5*4581 

-1.8156 

0. 

1.0098 

0 . 

-0.020521 

0. 

-0.23377 

0 . 

0.5383*4 


Table 15 

Exponential Coefficients Resulting from least 
Squares Fit, Stage it— 5* Based Upon 62 Location Years 



Const. Only 

Lin Only 

Lin T„ Only 
w 

Lin Only 

H 1 

-2.745776 

-U. 628175 

-3.243995 

-2.913867 

h 2 (xio' 2 ) 

0 . 

6. 3364 

0 . 

0 . 

h 3 (xkT 2 ) 

0 . 

0 . 

3-3732 

0 . 

Hj^XlCf 2 ) 

0 . 

0 . 

0 . 

5.9690 

H 5 (xicf 2 ) 

0 . 

0 . 

0 . 

0 . 


T x* D 1 & T A 

V d l & 'W 

Lin & Sq v/o P 

Lin & Sq v/ P 

H 1 

-9.155085 

-3.615012 

-17.68370 

-14.08316 

H 2 (X10* 2 ) 

23.795 

0 . 

53.654 

28.004 

H 3 (X10 -2 ) 

0 . 

5.1311 

33-366 

30.086 

H^CXIO -2 } 

156.75 

13.258 

261.05 

323-03 

H 5 (X10~ 2 ) 

0 . 

0 . 

0 . 

-2.4198 

h 6 (xio -3 ) 

0 . 

0 . 

-8.0340 

-3.9726 

H^XIO -3 ) 

0 . 

0 . 

-11-335 

-9.9229 

Hg(X10~ 3 ) 

0 . 

0 . 

-454.41 

-591-87 

H q (X10 -3 ) 

0 . 

0 . 

0 . 

1.0829 


-6.0586 

0 . 

0 . 

0 . 

H U (X10- ) 

0 . 

-0.62985 

0 . 

0 . 

H 12 (xio d ) 

0 . 

0 . 

0 . 

0 . 

h 13 (xio‘ 2 ) 

0 . 

0 . 

0 . 

0. 

H llt (X10 -2 ) 

0 . 

0 . 

0 . 

0 . 

h 15 (xio' 2 ) 

0 . 

0 . 

0 . 

0 . 


Lin P r Only 

Lm T x , T m , 

Lin T x’ V d l* 

P 

r 

-2.725344 

-4.725186 

-4.496764 


0 . 

5-4513 

4.4536 


0 . 

2.6297 

3-5240 


0 . 

-1.0206 

-1.5875 


-0.58432 

0 . 

-1.3646 


Lin & Int v/o P y 

Lin & Int v/ P 

T 

All v/o P_ 

XT' 

All v/ P r 

-1.652579 

8.809258 

10.10662 

-7-448562 

-7-9340 

-25.578 

-25-278 

46.310 

18.101 

-52.169 

13.198 

-38.699 

-158.898 

-455.43 

-796.65 

-119-32 

0 . 

96.168 

0 . 

101.24 

0 . 

0. 

-18.315 

-10.860 

0. 

0. 

-15-801 

3.8181 

0 . 

0. 

-806.70 

-273.37 

0 . 

0. 

0 . 

7.3280 

6.4918 

9.2975 

43-473 

6.3837 

-2.0119 

19.130 

-1.3173 

8.8611 

0 . 

-23.308 

0 . 

-15-033 

-0.34727 

0.20559 

1.2506 

0.48109 

0 . 

-0.30876 

0. 

0.86107 

0 . 

-1.4927 

0. 

-2.7542 
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*Uses form of Equation 5. All other stages use form of Equation 
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TEST PROCEDURES 

In attempting a preliminary assessment of results in the preceding phase, 
a simple test program was written which applied the parameters generated to 
daily meteorological data values. Erratic behavior was noted in using this 
procedure for the more highly nonlinear functions. This was believed due to 
the fact that, due to the nature of the phenological reports, only meteoro- 
logical values averaged over the phenological period were feasible for inclu- 
sion in the least squares parameter estimation. When parameters estimated on 
this basis were applied to daily meteorological values, many of these weather 
variables were out. ue the range used in parameter estimation, and great 
fluctuations were often the result for the nonlinear functions . 

When simpler functional fonns were utilized, these large fluctuations 
were not noted. Wien linear functional forms were utilized, reasonable fits 
were produced. Figures 5 through 11 depict the "tracking" of linear functional 
forms driven by daily meteorological values for typical functional forms and 
location-years . 

In an attempt to alleviate these difficulties, test programs were gen- 
erated using averaged meteorological values to drive the model on a daily 
basis. These were of two types. 

Interval Tests 

In these tests, various averaging techniques were applied to the inde- 
pendent meteorological variables used to drive the model, and statistics 
computed independently within each stage Interval, The data used in these 
tests was of necessity limited to those location-years for which reports were 
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available for both the beginning and ending date for the stage. This had the 
effect of limiting the size of the' data sets, particularly for stage intervals 
0-1 and 1-2, since only a few states adequately report emergence. A summary 
of the number of location-years included for each stage is shown: 

Stage Location-Years 

Plant-Eraerge (0-l) 18 

Emerge- Joint (1-2) 16 

Joint-Head (2-3) 33 

Head-Soft Dough (3-4) 30 

Soft Dough-Ripe (4-5) 25 

Specific location-years included in this interval test data set are summarized 

in Table 17. 

Running Averages . In this approach, the meteorological variables used 
to drive the model for any one day would be given by 


x r 



i=,1-n+I 

n 


Eq. 30 


in which x represents any of the independent variables T„, T„, D T , P . The 
value n is the number of days over which running averages are to b- computed, 
and may be varied within the program from 1 to 9, J is the Julian day number 
for which the growth rate is to be computed, and x' represents the averaged 
variable value to be used in the equation to compute growth rate for that day. 


HO 


Table 17 

Summary of Location Years Used 
for Interval Test Procrams 


Starve 


0-1 


1-2 


2-3 


State 

Crop Year 

C.R.D. 

Colorado 

1973 

6 


1975 

7 

Idaho 

1974 

9 

Oklahoma 

1964 

1 


1967 

1,2, 3, 4, 5, 6 


1969 

9 


1970 

9 


1971 

9 

Texas 

1975 

1S,1N,2N,8N, 

Colorado 

1973 

6 


■ 1974 

7 

Idaho 

1974 

9 

Oklahoma 

1965 

1 


1967 

1,2, 3, 4, 5, 6 


1969 

9 


1970 

9 


1971 

9 

Texas 

1975 

1S,2N,5 

Colorado 

1973 

6 


1974 

7 

Idaho 

1974 

1,9 

Kansas 

1964 

5 


1967 

1,9 


1969 

1 


1973 

5 


1974 

9 


1975 

7 

Montana 

1973 

1,2,3 

North Dakota 

1973 

1,4, 5, 6, 7 ,9 
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Table IT 
(Continued) 


Stap;e 

2-3 


3-4 


4-5 


State 

Cron Year 

C.R.D. 

Oklahoma 

1901 

1 ,2,3 ,4 ,5 


1 909 

9 


1970 

9 


1971 

9 

Texas 

1975 

1S,2N,5 

Colorado 

1974 

7 

Idaho 

1974 

1,9 

Kansas 

1964 

5 


1967 

1,9 


1969 

1 


1973 

5 


1974 

9 


1975 

7 

Missouri 

1971 

9 


• 1973 

9 


Montana 

1973 

1,2,3 

Oklahoma 

1967 

1,2, 3, 4, 5 


1969 

9 


1970 

9 


1971 

9 

Texas 

1975 

1S,2S,2N,5 

Idaho 

1974 

1,9 

Kansas 

1964 

5 


1967 

1,9 


1969 

1 


1974 

9 


1975 

7 

Missouri 

1973 

9 

Montana 

1973 

1,2,3 

Oklahoma 

1967 

1,2, 3, 4, 5 


1969 

9 


1970 

9 


1971 

9 

Texas 

1975 

IS ,2S ,2N ,! 
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Thus 


R « 


f (H,X') 


Eq. 31 


in which the function f may represent the polynomial, exponential, or 
Robertson equations discussed earlier. Within any interval the stage would 
be computed as 


J 

S„ * (K - 1) + i E i Eq. 32 

K 1 ' 

in which K represents the interval, represents the reported date at which 
the interval began, and J represents the Julian day number at which the stage 
is computed. 

Accumulative Averages . In this approach, meteorological values used to 
drive the model for any one day were computed as accumulative averages from 
the beginning day for the stage up to and including the day itself. Therefore, 
a typical averaged weather variable would be computed for a date j , as 


J 

I*i 

i=J G 

x’ * 2 Eq. 33 


in which x' represents the averaged weather variable used to compute growth rate 
for Julian date J, J o represents the Julian date at the beginning of the 
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phenoloGical Interval, and n represents the number of days over which averaging 
was done and would be computed by n * J - + 1. 

The rate for each day would then be computed by 

R * f (H,X’ ) Eq. 3U 

and Intermediate stage values within an interval would be computed by 

S R - (K - 1) + nR, Eq. 35 

Special considerations are necessary in applying this technique over 
interval 1-2, which includes a winter dormancy period. If the method were 
applied directly, then high temperatures and daylengths carried from the fall 
period would have a tendency to mask the onset of winter and subsequent 
dormancy. A similar problem would occur in the spring greenup period, when 
low values included from the winter would -mask the onset of spring and sub- 
sequent greenup. Thus, the winter dormancy and spring greenup criteria 
discussed earlier was used. The averaging algorithm was re-started when 
winter dormancy was indicated (10 days in which average 'ueraperature was less 
than k.69°C.), and again at spring greenup (10 days greater than 6.lh°C.), 

Test Results . Both of the averaging techniques were applied to indepen- 
dent data on a daily basis. Results from the running average used with poly- 
nomial functions (TEST CCll) and used with exponential functions (TEST CC12) 
are summarized in Tables 18 and 19- Results from application of the running 
average algorithm applied to the Robertson model are shown in Table 22. 

Results from applying the accuinulat ve average algorithm to polynomial 
(TEST CC21), exponential (TEST CC22), and Robertson functions are summarized 
in Tables 20, 21, and 22 respectively. 


A 


81 31^1 



Table 19 

Summary of Statistics from. Interval Test Program 
Using 5-Day Running Averages and Exponential Functions (TEST CC12) 



Stg. 

0-1 (18 L.Y. ) 

Stg. 1- 

-2 (16 L.Y.) 

Stg. 

2-3 (33 L.Y.) 

Stg. 3-4 (30 L.Y.) 

Stg. 

4-5 (25 L.Y.) 

FUUCTIOM 

BIAS 

VAR. 

RMSE 

BIAS 

VAF.. 

RMSE 

BIAS 

YAH. 

RISE 

BIAS 

VAR. 

RMSE 

BIAS 

VAR. 

RMSE 

Const . Only 

+2.9 

45.8 

6.8 

-27.4 

1403-0 

37-5 

+2-5 

■ 63-8 

8.0 

+2.1 

38.7 

6.2 

+9.0 

138.0 

11-7 

Lin V Only 

+2-9 

56.1 

7.5 

+2.2 

229-9 

15-1 

+1.9 

55-4 

7-4 

+1.7 

34.4 

5-9 

+8.0 

119-9 

11.0 

Lin Ty Only 

+2.2 

41.9 

6.5 

-3.6 

228.8 

15-1 

+1.6 

55-3 

7-4 

+1.2 

31-8 

5-6 

+8.2 

124.1 

11.1 

Lin C L Only 

+2.9 

1*5-8 

6.8 

-31.4 

1491-9 

38.6 

+1-9 

62.9 

7-9 

+2.6 

54.3 

7-4 

+8.7 

128.8 

11.3 

Lin P Only 

T 

+2.4 

1*3.7 

6.6 

-9-2 

970.3 

31-1 

+2.5 

63.5 

8.0 

+1.4 

35-0 

5-9 

+9-0 

138.2 

11.8 

Ian T x , T m> D l 

+2.3 

64.2 

8.0 

+3-3 

218.0 

14.7 

+1-5 

58.6 

7-6 

+2.7 

55-3 

7-4 

+7.6 

115.8 

10.8 

Lin T.^, T,_j, D^, V ^ 

+2.6 

69-2 

8.3 

+4.2 

249-8 

15-8 

+1-5 

55-2 

7.4 

+2.3 

46.5 

6.8 

- +7-7 

119-0 

10.9 

V °L & .** 

+2.9 

56.1 

7-5 

+30.6 

26.17.1 

51.2 

+1.4 

51.2 

7-2 

+1.0 

40.2 

6.3 

+8.7 

144.3 

12.0 

t m * d l & Vl 

+2.2 

1*1.9 

6.5 

+6.6 

1320.6 

36.3 

+1.0 

55-5 

7.4 

+1.1 

37-6 

6.1 

+7-9 

119-4 

10.9 

Lin h. Sq v/o P 

r 

+2.8 

63.1* 

8.0 

-8.0 

208.2 

14.4 

+6.1 

121.4 

11.0 

+5-5 

84.9 

9-2 

+4.8 

75-0 

8-7 

Lin & Sq v/ P 

r 

+2.? 

67-1* 

8.2 

-10-9 

281.1 

16.8 

+4.0 

77-9 

8.8 

+5-4 

79-2 

8.9 

+5-4 

88.4 

9-4 

Lin & Int w/o P 

r 

+2.6 

54.2 

7-4 

+6.4 

344.9 

18.6 

+0.2 

125-8 

11.2 

+0.5 

4l.4 

6.4 

+8.4 

139.6 

11.8 

Lin & Int v/ P 

r 

+1.8 

69.2 

8.3 

+116.0 

17502.3 

132.3 

+6.4 

118.3 

10-9 

-3-3 

41.6 

6.4 

-3.5 

117-3 

10.8 

All v/o P r 

+3.6 

67.1* 

8.2 

*+1.2 

298.3 

17-3 

+8.2 

170.8 

13.1 

+3.5 

74.0 

8.6 

+4.4 

1G0.2 

10.0 

All v/ P 
r 

+1.3 

86.8 

9-3 

125-8 

18193.7 

134.9 

-63.3 

11405-7 

106.8 

+11.4 

156.9 

12.5 

-35-7 

18349.4 

135-5 


TeCole 20 

Summary of Statistics from Interval 
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Stg. 0-1 (18 L.Y. ) 


Stg. 4-5 (25 L.Y.) 


Table 21 

Summary of Statistics from Interval Test Program 
Using Accumulative Averages and Exponential Functions (TEST CC22) 


Stg. 1-2 (16 L.Y.) 


FUNCTION 

BIAS 

VAB. 

RMSE 

BIAS 

VAR 

RMSE 

Const . Only 

+2.8 

42-9 

6.5 

-27-4 

1403.3 

37-5 

Lin T^ Only 

+2.5 

53.7 

7.3 

+0-9 

617.4 

24.8 

Lin Tj Only 

+2.1 

39-4 

6.3 

-9-4 

683-2 

26.1 

Lin D Only 

i_j 

+2.8 

42.9 

6.5 

-42.6 

2400.0 

48.9 

Lin P r Only 

+2.2 

42.6 

6.5 

-16.3 

1271.4 

35.7 

Lin T x , T m , D l 

+1.8 

66.4 

8.2 

+3.8 

620.4 

24.9 

Lin T x , T m . D l , P r 

+1.9 

69.4 

8.3 

+5.8 

699.2 

26.4 

t y> d l & t x d l 

+2.6 

54.8 

7-4 

+36.9 

4587.3 

67.7 

T M’ d l & t m d l 

+2.1 

39-4 

6.3 

+0. 4 

2762.9 

52.6 

Lin A Sq w/o P 

r 

+2.2 

60.8 

7-8 

-0-9 

389-1 

19-7 

Lin A Sq y/ 

+2.3 

65.4 

8.1 

-1.6 

350.2 

18.7 

Lin A Int w/o P_ 
r 

+1.8 

66.4 

8.2 

+13-1 

1138.2 

33-7 

Lin A Int w/ P_ 
r 

+2.1 

67.3 

8.2 

+28.3 

3319-1 

57-6 

All w/o P_ 
r 

+2.8 

60.9 

7-8 

+1.0 

615.0 

24.8 

All w/ P r 

+4.3 

76.6 

8.8 

+16.3 

2008.7 

44.8 


Stg. 2-3 (33 L.Y.) Stg. 3-4 (30 L.Y.) 


BIAS 

VAR. 

RMSE 

BIAS 

VAR. 

RMSE 

BIAS 

VAR. 

R1-5SE 

+2.6 

62.8 

7-9 

+2.1 

38.4 

6.2 

+9-0 

138.9 

11.8 

+1.7 

52.3 

7-2 

+1.6 

34.3 

5.9 

+7.9 

125.9 

11.2 

+1.2 

54.1 

7-4 

+1.1 

30.6 

5-5 

+8.0 

122.7 

11.1 

+1.6 

58.1 

7-6 

+2.6 

56.8 

7-5 

+8.7 

129.8 

11.4 

+ 2.7 

77.4 

' 8.8 

+1.1 

33-9 

5-8 

+9-0 

138.4 

11.8 

+1.1 

55.1 

7.4 

+ 2.7 

58.3 

7-6 

+7-5 

119.3 

10.9 

+1.0 

56.5 

7-5 

+2.3 

47.3 

6.9 

+7:6 

123.2 

11.1 

+0-5 

47.0 

6.9 

+1.0 

41.3 

6.4 

+8.7 

154.8 

12.4 

+0.5 

50.3 

7-1 

+0-9 

36.6 

6.1 

+8.0 

119-8 

10 .9 

+5-8 

175.7 

13.3 

+6.0 

95-9 

9-8 

+4.2 

83.8 

9-2 

+2.8 

117.7 

10.8 

+5-7 

87.4 

9-3 

+4-7 

92.7 

9-6 

- 0.1 

53-9 

7-3 

-1.4 

39-6 

6.3 

+7-5 

118.4 

10.9 

+2.5 

100.0 

10.0 

-1.6 

23-6 

4-9 

+7.8 

153.0 

12.3 

+7.8 

233.4 

15-3 

+ 4.9 

79-9 

8.9 

+2.8 

112.5 

10.6 

+4.7 

191.3 

13.8 

+ 8.9 

125-6 

11.2 

+7.0 

125-9 

11.2 


vn 

vn 


5 6 


Table 22 

Summary of Statistics from Interval 
Test Programs for Robertson Model 


Using 5-Day Running Averages 


Stage 

Bias 

Variance 

R.M.S .E. 

0-1 

+2.4 

66.8 

8.2 

1-2 

+1.1 

l4l. 3 

11.9 

2-3 

+1.4 

60. 8 

7.8 

3-4 

+2.3 

49-0 

7.0 

4-5 

+7.6 

115.2 

10.7 


Using Accumulative 

Averages 


Stage 

Bias 

Variance 

RaMaSaHj* 

0-1 

+2.0 

64.5 

8.0 

1-2 

-1.4 

551.0 

23.5 

2-3 

+1.2 

61.8 

7-9 

3-4 

+2.3 

50.0 

7.1 

4-5 

+7.2 

120.9 

11.0 


57 

Crop Year Teats 


In applying the above testing- procedures, data were limited within each 
stage to location-years for which both beginning and ending dates were reported 
for the stage. In order to overcome this limiting aspect as to location-years, 
a crop year test progi’am was written, in which the equations generated from 
least squares were applied throughout the entire winter wheat crop year, 

In this way, cumulative error throughout the growing season could be 
estimated. Further location-years for which phenological dates were not known 
in every interval could be included. For those intervals not having known 
phenological dates, the unknown error when a stage was reached was carried for- 
ward, and included within the cumulati /e error at the next stage at which a 
phenological date was reported. This program utilized the running average 
technique. 

A data set, consisting of 28 location-years of complete crop year meteoro- 
logical values was prepared for this program. The specific location-years 
included, and the stages reported for each, are shown in Table 23. 

Results of these tests are shown in Tables 24 and 25. The Robertson 
and a "combination" model are included in Table 24, In general, the crop 
year tests resulted in errors less than would be expected from the interval 
tests. The Robertson model, in particular, resulted in .relatively small 
biases and variances. 

A virtually infinite number of combinations, utilizing different func- 
tional forms may be investigated using this program. The results from one such 
combination, selected based upon bias and variance values from interval tests, 
are shown at the bottom of Table 24. This combination consists of the 
following functions; 
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Table 23 

Location Years Included in 
Crop Year Test Data Set 


State 

Colorado 

Idaho 

Kansas 

Missouri 

Montana 

North Dakota 


Oklahoma 


Texas 


C.R.D. 

Cron Year 

Stages Renorted 

6 

1973 

1,2, 3, 4, 5 

1 

1974 

2, 3, 4, 5 

9 

1974 

1,2, 3, 4, 5 

1 

1967,1969 

2, 3, 4, 5 

7 

1975 

2, 3, 4, 5 

9 

1967,1974 

2, 3, 4, 5 

9 

1971,1973 

3,4,5 

9 

1974 

3,5 

1 

1973 

2, 3, 4, 5 

2 

1973 

2, 3, 4, 5 

3 

1973 

2, 3, 4, 5 

1 

1973 

2,3,5 


1973 

2,3,5 

6 

1973 

2,3,5 

7 

1973 

2,3,5 

9 

1973 

2,3,5 

1 

1967 

. 1,2, 3, 4, 5 

3 

1967 . 

1,2, 3, 4, 5 

5 

1967 

1,2, 3, 4, 5 

9 

1969,1970,1971 

1,2, 3, 4, 5 

2N 

1975 

1,2, 3, 4, 5 

2S 

1975 

1,2, 3, 4, 5 

5 

1975 

1,2, 3, 4, 5 


Table 24 

Summary of Statistics from Yearly Test Program Using 
5-Bay Running Averages and Polynomial Functions {YEAR TEST l) (28 L.Y.) 

Stage 0-1 Stage 1-2 Stage 2-3 Stage 3-4 Stage 4-5 


FUNCTION 

BIAS 

VAR. 

RMSE 

BIAS 

VAR. 

RMSE 

BIAS 

VAR. 

RMSE 

BIAS 

VAR. 

RMSE 

BIAS 

VAR. 

RMSE 

Const. Only 

+1-9 

22.3 

U.7 

-2.1* 

11*60.7 

38.2 

-5.1* 

1308.0 

36.2 

-9-3 

897.3 

30.0 

+15-7 

862.0 

29-4 

Lin Only 

+1.3 

30.7 

5-5 

-8.5 

651.7 

25-5 

-9.8 

615-2 

21*. 8 

-2.5 

1*23-9 

20.6 

+1.5 

463-5 

21-5 

Lin T.j Only 

+1.4 

19.1* 

!*.l* 

— 1* . 3 

316.3 

17-8 

-3.6 

301.8 

17.1* 

+0.7 

216.7 

14-7 

+5-0 

372*8 

19-3 

Lin D t Only 

-U 

+1.9 

22.3 

U.7 

-14.3 

959-0 

31.0 

-l6.2 

90 U .1 

30.1 

-19.3 

761.6 

27-6 

+0.4 

219-5 

14.8 

Lin P r Only 

+1.1 

19-5 

!*.l* 

+5.3 

871 *. 0 

29-6 

+10.5 

712.7 

26.7 

■ +6.1 

578.6 

24.1 

+ 18.8 

789-6 

28.1 

Lin T X , T m , D l 

+0.6 

29.5 

5.1* 

-6.5 

1 * 91 *. 0 

22.2 

-8.3 

387.0 

19-7 

-1.3 

280.8 

16.8 

0.0 

424.2 

20.6 

Lin T x , T m , D l> P r 

+0.6 

32.1 

5.7 

-6.1* 

5l*8.8 

23.1* 

-8.7 

1*31-9 

20.8 

-1.8 

305.3 

17.5 

+0.6 

443-0 

21.0 

T x’ d l & ¥l 

+1.3 

30-7 

5.5 

-12.9 

1208.9 

31*. 8 

-11.1* 

670.0 

25.9 

-2.9 

11*8.8 

12.2 

+4.3 

235-4 

15-3 

V D L & Vl 

+1.1* 

19- 1 * 

!*.!* 

-6.3 

591-1* 

2l*.3 

-9-9 

381*. 2 

19.6 

-3.6 

232.1* 

15.2 

+7.1 

300.7 

17-3 

Lin & Sq w/o P r 

+1.5 

23-6 

k.9 

-20.7 

11*31.1 

37-8 

-l4.8 

1007-1* 

31.7 

+ 0 . 1 * 

543.1 

23.3 

+9-5 

267.1 

16.3 

Lin & Sq v/ P^ 

+1.2 

25.5 

5.1 

-28.1 

1921.3 

1*3-8 

-17.1* 

871.2 

29.5 

-3.2 

252.2 

15.9 

+6.8 

248.9 

15-8 

Lin & Int v/o P 

r 

+0.5 

29.6 

5-1* 

-14.8 

1153.7 

3U.« 

-13.1 

1109.1* 

33-3 

-0.3 

252.6 

15-9 

+8.5 

444.1 

21.1 

Lin & Int v/ P 

r 

+0.8 

31-1* 

5.6 

-7-5 

1736.3 

1*1.7 

+5.1* 

987.1* 

31.1* 

-2.5 

452.1 

21.3 

+11.4 

858.5 

29.3 

All w/c P r 

+2.1 

20.6 

1*. 5 

-9-7 

935.5 

30.6 

+9-9 

1178.9 

34.3 

+ 26. 8 

1725.0 

41.5 

+44.4 

3571-6 

59-8 

All v/ P 

r 

+3.2 

27-9 

5-3 

-10.3 

1069.3 

32.7 

-11.8 

889.3 

29.8 

-5-2 

362.3 

19-0 

+ 4.9 

389-9 

19-7 

Robertson 

+1.5 

26.3 

.5.1 

-2.0 

258.5 

16.0 

« 

-3-3 

200.5 

14.2 

-0.2 

205.6 

14.3 

+1.8 

305.5 

17-5 

Combination 

+1.1* 

19-1* 

l*.i* 

-3-1* 

262.8 

16.2 

-5.5 

186.6 

13.7 

-2.7 

125-7 

11.2 

-2.9 

197-8 

14.1 


Table 25 

Summary of Statistics from Yearly Test Program Using 
5-Day Sunning Averages and Exponential Functions (YEAR TEST 2) (28 L.Y. ) 

Stage 0-1 Stage 1-2 Stage 2-3 Stage 3-4 Stage 4-5 


FUNCTION 

BIAS 

VAR. 

RMSE 

BIAS 

VAR. 

RMSE 

BIAS 

VAR. 

RMSE 

BIAS 

VAR. 

RMSE 

BIAS 

VAR. 

RMSE 

Const. Only 

+1.9 

22.3 

4-7 

-1.4 

1456.8 

38.2 

-5.4 

1308.0 

36.2 

-9-3 

897-3 

30.0 

+14.7 

831-7 

28.8 

Lin Tjj. Only 

+1.5 

28.9 

5.4 

-1.6 

348.5 

18.7 

-3-8 

334.4 

18.3 

-.1 

320.1 

17-9 

+2.7 

396.9 

19-9 

Lin Tjj Only 

+1.4 

19.4 

4.4 

-1.2 

199-1 

l4.1 

-1.0 

217-9 

l4.8 

-0.8 

198.1 

14.1 

+5-6 

294,3 

17.2 

Lin D_ Only 

±J 

+1.9 

22.3 

4.7 

-9.3 

1066.3 

32.7 

-11.9 

983.8 

31.4 

-15-9 

747-4 

27-3 

+6.2 

336.3 

18.3 

Lin P„ Only 

Y 

+1.3 

18.0 

4.2 

+13.8 

1098.6 

33.1 

+0.2 

4224.0 

65 .O 

+3.4 

2830.0 

53-2 

+36-5 

3049.1 

55-2 

Lin T x> Tjj, D^ 

+0.8 

27.2 

5.2 

-1.8 

356.7 

18.9 

-3.4 

263.8 

16.2 

+0.9 

321.6 

17-9 

+3-7 

397-0 

19-9 

Lin T x , Tjj, D^, P r 

+1.0 

27-9 

5-3 

-1.2 

429-0 

20.7 

-3-5 

317-7 

17.8 

+0.9 

372.6 

19-3 

+5-2 

415-6 

20.4 

V d l * Vl 

+1.5 

28.9 

5-4 

-13.4 

3457.1 

58.8 

-10.1 

877-3 

29-6 

- 0.6 

698.8 

26.4 

+8.0 

639-5 

25-3 

t m’ d l * t m°l 

+1.4 

19-4 

4.4 

-22.2 

4410.3 

66.4 

-5-2 

1198.6 

34.6 

+ 7.4 

1024.5 

32.0 

+13-6 

1079-7 

32*. 9 

Lin & Sq w/o P r 

+1.5 

24.2 

4.9 

-14.7 

1057-2 

32.5 

-3-3 

452.9 

21.3 

+5-9 

474.5 

21-7 

+4.7 

341.2 

18.5 

Lin 8t Sq w 1 P^ 

+1-5 

24.7 

5-0 

-9-8 

532.7 

23-1 

-7.6 

355-3 

18.8 

- 0.5 

332-6 

18.2 

+2.7 

281.6 

16.8 

Lin & Int w/o P 

r 

+0.8 

27-2 

5-2 

-10.1 

1127.8 

33.6 

-9.9 

976. ID 

31.2 

-1-9 

427.1 

20.7 

+ 4.7 

656.2 

25.6 

Lin & Int w/ P^ 

+1.1 

27-8 

5.3 

-38.0 

10075-8 

100.4 

-32.0 

8536.4 

92.4 

+ 2.9 

2133-4 

46.2 

+5-7 

3341-7 

57-8 

All v/o P r 

+•1.9 

21.0 

4.6 

-16.4 

1925.1 

43-9 

+3.8 

1038.7 

32.2 

+14.1 

1291.1 

35-9 

+25-5 

2100.3 

45.8 

All w/ P r 

+3.5 

30.4 

5-5 

-19-7 

5600.0 

74.8 

-7-9 

2083.5 

45.6 

+4.9 

1913.9 

43-7 

+13.7 

1872.3 

43.3 
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Stage 

Function 

Terms in Fc nation 

0-1 

Polynomial 

Linear T^ Only 

1-2 

Robertson 

Triquat ratio 

2-3 

Exponential 

T x , D l and T. J) l 

3-4 

Polynomial 

Linear T^ Only 

4-5 

Exponential 

Linear and Oq, w/o P . 

r 


This particular combination worked well, resulting in low variances, relative 
to the other cases tested, although biases exceeded slightly those from th'* 
Robertson model. 

Figure 12 shows, in graphical form, the general behavior of the bias 
values resulting from these crop year tests, for typical functions, and 
Figure 13 shows a similar plot ‘for variance values. 


VARIANCE PROPAGATION 

In order to assess the acceptability of the variances resulting from the 
testing procedures, the technique of variance propagation may be used as an 
independent check. The generalized variance propagation relation provides a 
powerful tool to determine the resulting variances of quantities computed from 
parameters and independent variables subject to statistical variation. 

Mathematical Model 

The basic relationship utilized in variance propagation is 



J 0 j* 

= SL 

c,n «,n n,c 


c,c 


Eq. 36 






(sjCbq) TS’WH 


Plot of Variance from Crop Year 
Test for Various Cases 



In this equation, represents a c x c variance-covariance matrix for the 
c x 1 vector of functions 


c,l c,l n,I * 


Eq. 37 


The n x 1 vector jr represents the quantities (parameters and independent 
variables) upon which computation of s_ is based, through the functional rela- 
tionships f, and Qy is the n x n variance-covariance matrix of these quantities. 
For the problem under consideration 


s 

1,1 


R 


f (JL JL 


) 


Eq. 38 



Eq. 39 


where n = p + q, H is the vector of parameters estimated from least squares, and 
X is the vector of independent weather related variables defined previously. 

Then Equation 3 6 may be applied to growth rate as 


„ = Jsz .. Q y... Lm 

q R 1,19 19,19 19,1 


Eq. 40 


where q^ is the variance of daily growth rate, and 
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Si 

0 

.A. . 

15,15 

157T 

19,19 

0 

Sc 


_ ¥7i5 

W J 


Eq, 1*1 


The key to application of the relationship is to obtain values for the parameter 
variance-covariance matrix C^. One of the attributes of the generalized least 
squares technique used is that propagation is carried out during parameter 
estimation, and reliable estimates of the parameter variance-covariance matrix 
result from the least squares estimation process (3,6). 

After estimation of the variance of a typical growth rate within a given 
stage, q R , other variances may be computed. The number of days within a stage 
may be approximated by 

A A 

n * l/R Eq. 42 

A 

in which n is the estimated number of days in that stage, and R is a typical 
growth rate within that stage given by 

R = f(H.X) Eq. 43 

A 

where X represent typical values of weather related variables within the stage, 
e.g. , average or median values. Using Equation 42 and applying the variance 
propagation relationship of Equation 36 results in 



Eq. 44 
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and a standard deviation of 

°n * ^ ' Eq. **5 

These represent the estimates of variance and standard deviation of the number 
of days within each stage, which may be compared to those values computed by 
the interval tests. To estimate the cumulative variance, in days, at the end 
of any stage j , the relationship 

i 

Nj = l Eq. U6 

J=1 

may be used. Applying the variance propagation relationship of Equation 36 
to this function results in 

J 

S * I < t ** 

i=l 1 


and 




‘Eq. H8 


Test Results 

The above relationships were programmed and run on the computer. Variance- 
covariance matrices for parameters were taken directly from the least squares 
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program. The interval test data set was used to compute typical weather-related 

variables, in this ca^e averages, as well as variances for and correlations 

between these variables. The results revealed a high correlation between the 

T„, T„, and D_ variables, and these correlations should not be neglected in 
XML 

further work. 

The results of these computations are, for selected cases, shown in 
Table 26. The figures shown are the variances and standard deviations pre- 
dicted by variance propagation for the cumulative mismatch in days, which may 
be compared directly to the results from the crop year test sequence. Figures 
l4, 15, l6, and 17 show these comparisons, in graphical form, for the polynomial, 
exponential, Robertson, and combination models, 


Wr- 




Table 26 

Cumulative Variances and Standard Deviations 
Predicted fey Variance Propagation for Selected Cases 


FUNCTION 

Stage 0-1 

Stage 1-2 

Stage 2-3 

Stage 3-4 

Stage 4-5 


Var 

(Days ) 

Std. Dev. 
(Days) 

Var ? 
(Days ) 

Std. Dev. 
(Days) 

Var „ 
(Days ) 

Std. Dev. 
(Days) 

Var_ 
(Days ) 

Std. Dev. 
(Days ) 

Var_ 
(Days ) 

Std. Dev. 
(Days) 

Polynomial 











Const . Only 

7-34 

2.71 

125.66 

11.21 

125.66 

11.21 

125.89 

11.22 

125.95 

11.22 

Lin T„ Only 
/* 

4.71 

2.17 

660.49 

25-70 

662.03 

25-73 

662.03 

25.73 

666.16 

25.81 

Lin Only 

3.67 

1-92 

712.89 

26.70 

716.63 

26.77 

741.47 

27.23 

742.02 

27 24 

V D L 4 T A 

4.71 

2.17 

449.86 

21.21 

501.31 

22.39 

501.31 

22.39 

503.55 

22.44 

V d . & Vl 

3-67 

1.92 

454.54 

21.32 

468.25 

21.64 

468.72 

21.65 

470.02 

21.68 

Lin T x , T m , D l 

3.87 

1.97 

506.25 

22.85 

523.49 

22.88 

523.95 

22.89 

529.92 

23-02 

Lin T^, T,j, D^, P f 

7-89 

2.81 

924.77 

30.41 

974.06 

31.21 

979.69 

31.30 

1000.46 

31.63 

Lin & Set v/o P^ 

8.05 

2.84 

338.56 

18.40 

349.69 

18.70 

352.31 

18.77 

355-70 

18.86 

Lin 4 Si s/ P 

r 

13.66 

3.69 

1730.56 

4l.6o 

2043. 04 

45.20 

2055.72 

45.34 

2058.44 

45-37 


c\ 

C3 


Table 2 6 
( Continued) 


FUNCTION 

Var 

Std. Dev. 

Var_ 

Std. Dev. 

Var„ 

Std. Dev. 

Var_ 

Std. Dev. 

Var_ 

Std. Dev. 

* 

(Days ) 

(Days) 

(Days ) 

(Days) 

(Days ) 

(Days) 

(Days ) 

(Days) 

(Days ) 

(Days) 

Exponential 

Const. Only 

7-35 

2.71 

126.11 

11.23 

126.34 

11.24 

126.56 

11.25 

126.70 

11.25 

Lin T^ Only 

5-30 

2.30 

1073.87 

32.77 

1075.18 

32.79 

1076.00 

32.79 

1079.78 

32.86 

Lin T„ Only' 

3.69 

1.92 

1086.36 

32.96 

1090.98 

33-03 

1091.64 

33.04 

1092.39 

33-05 

T x* d l 4 T -A 

5-30 

2.30 

3237.61 

56.90 

3242.16 

56.94 

324.44 

56.96 

3249-01 

57-00 

V d l 4 T A 

3.69 

1.92 

3099-15 

55.67 

3103-60 

55.71 

3209.22 

56.65 

3210.36 

56.66 

Lin T x , T j5> B. 

3-68 

1.92 

921.12 

30.35 

922.94 

30.38 

923.02 

30.38 

929.64 

30.49 

Lm T’ x , T m , D l , P r 

5-86 

2.42 

1357.19 

36.84 

1395.02 

37.35 

1400.26 

37-42 

1405-50 

37.49 

Lin !t Sq. v/o P 

r 

8.58 

2.93 

1042. 6U 

32.29 

1047.17 

32.36 

1053.65 

32.46 

1055.60 

32.46 

Lin & Sq v/ P^ 

17. 1*1* 

4.18 

1031.05 

32.11 

1060.15 

32.56 

1064.06 

32.62 

1068.64 

32.69 

Bobertson 

35-90 

5. 59 

1046.52 

32.35 

1048.46 

32.38 

1050.41 

32.41 

1069.94 

32-71 

Combination 

3.67 

1.92 

992.25 

31.50 

997-30 

31.58 

1022.08 

31-97 

1025.28 

32.02 



R.M.S.E. from Year Test Compared to those Predicted 
from Variance Propagation for Typical Polynomial 




R.M.S.E. from Year Test Compared to those Predicted 
from Variance Propagation for Typical Exponential 








Figure 17 

R.M.S.E. from Year Test Compared to Those Predicted 
from Variance Propagation for Combination Model 




SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


Summary 

Winter wheat phenological information was gathered for crop reporting 
districts located in the Great Plains and Rocky Mountain regions of the 
United States. Corresponding meteorological data were obtained from the 
National Climatological Center for these location-years, for values of daily 
maximum temperature, minimum temperature, and precipitation. Daylength 
values were computed. 

These data, after editing and reduction, were utilized in least squares 
fitting programs using the method of simultaneous adjustment of observations 
and parameters (generalized least squares). Parameters were estimated using 
these techniques for various forms of polynomial, exponential, and the 
Robertson "triquadratic” functions. The functions resulting from these fits 
were then tested in a predictive mode by applying them to data sets withheld 
from the least squares procedures. Variance propagation studies were con- 
ducted in order to predict expected variances of growth rates computed from 
the parameters estimated. 

Conclusions 

The investigations conducted support the following general conclu- 
sions: 

1. ) The use of phenological reports from crop reporting 

districts was troublesome. Although use of such informa- 
tion requires no separate data gathering, as it is accom- 
plished within the U.S.D.A. reporting service, the data 
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have several disadvantages, Few states report all pheno- 
logical stages, and many report other than the ’’standard" 
stages (plant, emerge, Joint, head, soft dough, ripe). No 
states, to this investigator’s knowledge, report informa- 
tion concerning dormancy or spring greenup. 

2. ) The data from crop reporting districts are ’’noisy", and in 

3ome instance contained outright mistakes. This caused re- 
latively high variances in parameter estimates, and, in turn, 
relatively large uncertainties in predicted phenological 
stages. 

3. ) Because only end points are available, at best, in the 

phenological reports, only average growth rates may be com- 
puted and utilized for parameter estimation within each stage 
interval for each location-year. This, in turn, necessitates 
that the corresponding environmental values be averaged 
spatially and temporally over the CRD and the stage interval, 

4. ) For the emerge to Joint interval special considerations are 

necessary. This interval includes the long winter dormancy 
period, and some criteria must be established in data pre- 
paration which reflects this fact. Later testing procedures 
indicated that the behavior of the resulting models were 
relatively insensitive to minor variations in the dormancy 
assumptions made. 

5. ) The generalized least squares modeling techniques, combined 

with numerical methods for linearization and function 


evaluation, appear to provide a powerful and flexible tool 
for parameter estimation within a variety of function types. 
While some difficulties were encountered in obtaining con- 
vergence for models with higher order terms, due to the 
approximations involved, the problems were relatively minor, 
and are certainly reconcilable. Fairly good parameter approx 

filiations are necessary in using the model. Further, the 
method lacks efficiency in testing of variance for inclu- 
sion of terms when compared with stepwise regression. 

However, the flexibility of the technique, including 
its ability to accommodate nonlinear models, and 
configurations within each model is a great advantage. 

The fact that the method recognizes and includes computa- 
tionally the statistical variability within all observed 
quantities, and minimizes with respect to all these, is 
felt to result in more realistic parameter estimates than 
those resulting from regression analysis; 

6 . ) In virtually no instance was the inclusion of precipitation 

within the model found to be significant. 

7. ' Tests of functions on independent data revealed that the 

more highly nonlinear functional forms produced very erratic 
results when driven with daily environmental variables. 

This was not unexpected, since parameter estimates were gen- 
erated using averaged data, and the daily values used for 
testing often lay outside the range of these averaged values, 
resulting in greatly magnified results for those points. 
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8. ) Attempts to alleviate the difficulties of 7.) above through 
the use of running and accumulative averages was only par- 
tially successful. 

9. ) Variance propagation techniques appear to be a useful device 
for predicting expected uncertainties. In this study, these 
techniques were used to verify that the uncertainties noted 
during testing were reasonable for the primary data set used 
to generate and test parameters. 

Recommendat ions 

The following recommendations are made based upon the investigations 
conducted. 

1. ) The best combination of functions tested was the following: 

R = 0.87870 - 0.00045898 T m 

R * 0.071784 (D l - 0.2796) • {(4.6869 x 10” l4 )(T x + 2.3876) 
(4.6618 x 10" 6 )(T x + 2.3876) 2 + (2.3943 x 10- 4 ) 

(T m + 2.3876) + (1.7055 x 10’ 5 )(T m + 2.3876) 2 } 

R ={exp -6.22319 + 0.12008 T x + 1.0337 (D l - 12) - 
0.038591 T x (D l - 12 )} 

R = 0.37980 + C ,0061145 t m 
R = exp{-17.6837 + 0.53654 T x + 0.333 66 T M + 

2.6l05(D L - 12) - 0.0080340 T x 2 - 0.011335 T^ 2 - 
0.4544l(D T - 12) 2 }. 

Lt 

2. ) Future attempts at parameter estimation should be based upon 

< better, less "noisy" data, particularly as concerns 


Stage 0-1 : 
Stage 1-2 : 

Stage 2-3 : 

Stage 3-4 : 
Stage 4-5: 
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phenolcgical reports, Current efforts at JSC to collect data 
bases formulated from intensive test sites and blind sites 
represent a desirable first step in this direction, 

3. ) The generalized least squares technique is recommended for 

parameter estimation of any sort, rather than the more standard 
linear regression techniques. The method can accomodate both 
linear and nonlinear model forms , and provides useful var- 
iance information concerning parameter estimates, 
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